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Exact transition probabilities in a 6-state Landau-Zener system with path interference 
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We identify a nontrivial multistate Landau-Zener model for which transition probabilities between any pair of 
diabatic states can be determined analytically and exactly. In the semiclassical picture, this model features the 
possibility of interference of different trajectories that connect the same initial and final states. Hence, transition 
probabilities are generally not described by the incoherent successive application of the Landau-Zener formula. 

We discuss reasons for integrability of this system and provide numerical tests of the suggested expression for 
the transition probability matrix. 
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FIG. 1. (Color online) Diabatic energy levels (straight black lines) and their pairwise couplings (parameters =b^i, z = 1,2,3) shown at 
corresponding level intersections. Red and Blue arrows illustrate trajectories that start at the level-1 and end either at the level-4 or at the 
level-3 with, respectively, destructive and constructive interference. 


The multistate version of the two-state Landau-Zener model is one of the most fundamental systems in nonstationary quantum 
mechanics Eia, which finds numerous applications in atomic and condensed matter physics El. It considers the evolution of 
N states described by the Schodinger equation with parameters that change linearly with time: 

Here, t/) is the state vector in a space of N states; A and B are constant Hermitian N x N matrices. One can always choose the, 
so-called, diabatic basis in which the matrix B is diagonal, and if any pair of its elements are degenerate then the corresponding 
off-diagonal element of the matrix A can be set to zero, that is 

Bij = Anm = 0 if [3fi = Prri’) 72 7^ 777/ G (l,...A/^). (2) 

Constant parameters Pi are called the slopes of diabatic levels', nonzero off-diagonal elements of the matrix A in the diabatic 
basis are called the coupling constants, and the diagonal elements of the Hamiltonian, Hu = Pd + where = An, are called 
the diabatic energies 01. 

The goal of the multistate Landau-Zener theory is to find the scattering N x N matrix S, whose element Snn' is the amplitude 
of the diabatic state n at t ^ +oc, given that at t ^ —oc the system was in the n'-th diabatic state. In most cases, only 
the related matrix P, Pn'^n = called the matrix of transition probabilities, is of interest. Generally, for N > 2, the 

analytical solution of the model 0 is unknown. Nevertheless, a number of exactly solvable cases with specific forms of matrices 
A and B have been derived EHI3. 

Several influential results in the multistate Landau-Zener theory have been originally proposed as conjectures based on intu¬ 
ition and numerical calculations. This includes the solution of the generalized bow-tie model lEl, Brundobler-Elser formula 
m, and the no-go rule Q. The subsequent search for the proofs (SI not only confirmed these conjectures but also advanced our 
understanding of solvable systems in nonstationary quantum mechanics beyond the multistate Landau-Zener model alia. 

In this article, we introduce one more system of the type 0/ which matrix of transition probabilities can be written explicitly. 
This model describes the system of = 6 interacting states with the Hamiltonian 
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where parameters ei, 62, gi, g2, Qs, P are arbitratry constants and t is time. 

It is convenient to illustrate parameters of this model on a diagram that plots diabatic energies as functions of time and marks 
the pairwise intersections of diabatic levels by corresponding coupling constants, as shown in Fig. For the Hamiltonian ([^, 
this diagram shows that the model describes two crossing bands of energy levels, with three levels in each band. 

When diabatic energy levels in each band are well separated from each other, i.e. ei, 62 ^ 5 ^ 1 , 5 ^ 2 , ^ 3 , one can justify the 
semiclassical approach to estimate transition amplitudes. During time intervals for which all levels are well separated, diabatic 
states are almost eigenstates of the Hamiltonian. The adiabatic theorem predicts then that transitions between such states are 
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exponentially suppressed. The adiabaticity condition is violated only for pairs of crossing diabatic states. After passing through 
such a crossing point, a system can remain on the initial level or turn to another one. Hence, one can visualize possible trajectories 
of the system as we illustrate in Fig.[2 The blue and the red color arrows indicate the semiclassically allowed trajectories that 
connect the initially populated state at the level-1 with the levels 3 and 4. Here we note that only trajectories that turn to the right 
or upward are allowed in Fig.l^in the semiclassical limit. 

The transition probability between a pair of states after a passage through a crossing point is described by the two-state 
Landau-Zener (LZ) formula. In our case, if a diabatic level crosses another one with a corresponding nonzero coupling Qj, 
j G (1, 2, 3), then with a probability 


(4) 

the system will remain in the same diabatic state after the level crossing. Respectively, the probability to turn to another level is 
9j = ^-Pj- 

By applying the LZ-formula to pairs of diabatic level crossings according to their appearance in a chronological order one can 
estimate the amplitude of a semiclassical trajectory. For example, consider the process in Fig.l^that starts at the level-2 and ends 
at the level-4. Semiclassically, there is only a single trajectory that connects these states. It starts at the level-2, passes through 
the intersection with the level -6 without changing the diabatic state, then it turns to the level-4 at the corresponding intersection 
and then goes through the intersection with the level-3 without changing the level. The amplitude of such a process is given by 

S42 = \/p 3 (l (5) 


where is the dynamic phase of a trajectory: 


4>d = J Wkit)t + ek(t))dt, (6) 

and is the “LZ-phase” that the system accumulates after transitions through the crossing points. Time-dependence of 
indexes in ([^ indicates that different parts of a trajectory generally correspond to system’s presence at levels with different 
values of the index k and consequently different values of parameters Pk G (0, P) and ek G (— 62 , 0, ei). Since there is only a 
single trajectory connecting levels 4 and 2, the phase factor does not influence the corresponding transition probability: 

-P2^4 =P3(1 -Pl)P2- (7) 

The situation becomes more complex when there are more than one trajectory connecting initial and final states. In our case, 
such trajectories start from either the level 1 or 6 and end at either the level 3 or 4. In an arbitrary nonintegrable model, relative 
dynamic phases of different trajectories change violently with changing parameters. However, the model with the Hamiltonian 
0 is tuned so that different trajectories connecting the same initial and final states acquire the same dynamic phase. Consider, 
e.g. transitions from the level-1 to the level-4. One can find that the blue-color trajectory in Fig. spends the same amount of 
time at the level-1 with the constant part of the diabatic energy equal to ei as the red trajectory spends at the level-4 with the 
same value of the constant part of the diabatic energy, etc. Another way to see that the dynamic phase is the same for both paths 
is to use the property that the dynamic phase difference between two trajectories is proportional to the area between them cni. 
In Fig. [2 the blue and the red trajectories enclose two plaquettes of equal size but different sign. Hence, despite the trajectory 
amplitude interference, the dynamic phase drops out of the final expression for the transition probabilities in the model®. 

The role of the LZ-phase in the independent crossing approximation was discussed by Demkov and Ostrovsky in ifT^ They 
showed that the effect of this phase can be trivially included by assuming that, at each pairwise crossing point, if a system 
remains at the initial diabatic level it does not gain an additional phase factor but, if it changes the diabatic state, the amplitude 
of the trajectory acquires an additional factor i. We found that this rule is not sufficient for our case because it was derived for 
real positive coupling constants. In the model in Fig.[2 all coupling constants appear with opposite sings at different intersection 
points. To account for this case, consider a simple Landau-Zener system: 

i-^a =/3ita+\g\e^‘^<>b, i-^b = ^2tb + \g\e~^‘^^a, ( 8 ) 

where we explicitly separated the absolute value and the phase of the coupling constant. By a change of variables, b ^ 

the system <0 transforms into the one with a real positive coupling: 

= I3ita + \g\b, i^^b = I32tb + \g\a. (9) 

This means that the transition amplitude from the diabatic state described by the amplitude a into the other diabatic state is 
different in models ^ and ([^ by a phase factor . In application to the multistate Landau-Zener models, this means that 
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FIG. 2. (Color online) (a-b) Typical dependence of adiabatic energies (eigenstates of the Hamiltonian matrix) on time t. The model in (a) 
corresponds to the Hamiltonian ^ and (b) corresponds to the Hamiltonian ( p^ . The case (a) features three exact level crossings at t = 0 
while these degeneracies are lifted in the case (b). (c-d) Transition probabilities from the level-1 as functions of the level separation for a = e, 
62 = 1.5e, gi = 0.25, g 2 = 0.3, gs = 0.35, /3 = 1. The case (c) corresponds to the model (a) and (d) corresponds to the model (b). 
(e-f) Transition probabilities from the level-2 as functions of the level separation for gi — g 2 — gs — 0.3, ei = 62 = e, /3 = 1. The case 
(e) corresponds to (a), and (f) corresponds to (b). All discrete points correspond to results of direct numerical simulations of the evolution 
from t = —1000 to t = 1000, and solid lines are theoretical predictions of Eq. Details of the numerical algorithm are discussed in the 
supplementary material for Ref. m. 
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each time a semiclassical trajectory switches between different diabatic levels, its amplitude acquires a phase factor , were 
(l)g is the phase of the coupling constant at a corresponding level intersection. Following this rule, we should assume for the 
model in Fig. that whenever a trajectory changes levels, it acquires a phase factor ±i for, respectively, positive and negative 
values of the couplings at intersections. 

Applying this rule to the transition from the level-1 to the level-4 we find two contributing amplitudes: 

= {-I)^i^\/{1-P2){l-P3){l-Pl)p2, 5'4i"® = -i^\/p2(l -Pl)(l -P3)(l -P2)- (10) 

The sum of those amplitudes is identically zero, i.e., trajectories interfere destructively so that the transition from the level-1 to 
the level-4 is forbidden: = 0. For the transition from the level-1 to the level-3 we find: 

'S'si'' = (-l)^*'‘y(l -P 2 )(l -P3)(l -Pl)(l -P 2 ), 5^1"® = -iVP 2 (l -Pl)(l -P3)P2- (11) 

Amplitudes ( pT) interfere constructively, so that Pi ^3 = (1 — P 3 )(l — pi). We are now in a position to summarize the prediction 
of the independent crossing approximation for the form of the transition probability matrix: 
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So far, our derivation of the transition probability matrix was purely semiclassical. It is expected to be predictive only for 
sufficiently large values of parameters ei and 62 . In what follows, we are going to claim that for the model with the Hamiltonian 
the matrix ( p^ is actually exact, i.e. it is valid for arbitrary values of all parameters. 

In order to understand what is special about the model 0 , it is instructive to compare this model with a similar one that has 
the Hamiltonian 
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(13) 


interference and cancellation of the relative dynamic phase equally applies to the model ( p^ . For example, transitions that 
correspond to a single allowed semiclassical path, such as from the level -2 to any other state, should have the same probabilities 
in both models. The only semiclassical difference is that now constructive and destructive interference correspond to different 
transitions. 

Nevertheless, models ^ and ( p^ show very different behavior beyond the semiclassical limit. The crucial difference between 
them becomes apparent if we plot eigenstates of their Hamiltonians (adiabatic energies) as functions of the parameter t, as shown 
in Figs.[^a,b). Figure [^a) corresponds to the model It demonstrates that at f = 0 the system experiences three pairs of 
exact crossings of adiabatic energies. In contrast. Figure [^b) shows that, in the model ( p^ , those pairwise degeneracies are 
lifted. Returning to the diabatic basis, we find that the moment t = 0 also corresponds to the crossings between diabatic levels: 
1 and 4, 2 and 5, 3 and 6 . Such crossing levels are not coupled to each other directly, however, higher order processes can and do 
induce transitions between them. Nevertheless, the symmetry of the model ^ is sufficient to completely decouple those pairs of 


states at f = 0 but in the model (13) mini-gaps open up. Such mini-gaps are the signatures of the breakdown of the semiclassical 
approximation. 

Let us compare now the transition probabilities in models Q and ( p^ as functions of the level separation, as we illustrate 


m Figs. I^c-f) for different initial conditions and different ratios of ei/e 2 . In the case of the model numerically obtained 
transition probabilities are found to be independent of the values of the parameters ei and 62 , and they perfectly agree with the 
theoretical prediction ( p^ . In contrast, in the model ( p^ , transition probabilities generally depend on ei, and 62 nontrivially and 
saturate only in the limit ei, 62 ^ 91 ^ 92 , 93 , as it is expected from the semiclassical approximation. 

Comparison of the models 0 and ( p3] ) does support the conjecture that the transition probability matrix is the exact 
nonperturbative solution of the multistate Landau-Zener problem with the Hamiltonian As an additional numerical proof, 
we provide Fig.j^that compares the theoretical prediction ( p^ with numerical simulations of the model ^ with different choices 
of the coupling constants. All our numerical tests have shown perfect agreement with Eq. ( p^ . 

In conclusion, we identified and explored the Landau-Zener-like model of two crossing bands, for which we determined the 
matrix of transition probabilities. This model shows a relatively complex behavior due to the possibility of semiclassical path 
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FIG. 3. (Color online) Numerical test of Eq. for initially populated (a) level-1 and (b) level-5. All discrete points correspond to results of 

direct numerical simulations of the evolution with the Hamiltonian from t — —1000 iot — 1000. Solid lines are theoretical predictions of 
Eq. |T^ . The choice of parameters is: ei = 0.25, 62 = 0.35, gi = 0.85^, ^2 = 1.0^, = 1.15^, /3 = 1. 


interference leading to either constructive or destructive interference. Our numerical simulations confirm that the transition 
probability matrix, which we derived in a semiclassical framework, is exact i.e. it describes an arbitrary choice of the model 
parameters. 

Our model was identified on the basis of its similarity with other known exactly solvable systems, such as the reducible 6-state 
model in and the generalized bow-tie model do). For such previously solved systems, the semiclassically derived and exact 
expressions for transition probabilities coincide. Two properties are also shared by our model with these systems: (i) the absence 
of the dynamic phase effect on transition probabilities in the semiclassical framework and (ii) the exact crossing of adiabatic 
energies at intersections of diabatic states without direct couplings. This finding suggests that our model and the previously 
solved systems lH [Tol are likely the instances of a more general class of solvable models that share the same properties. It 
should be interesting to identify new candidate systems of the type of Eq. Q with properties (i-ii) and test them for integrability. 
Further progress in this direction will be also achieved by deriving a mathematically rigorous proof of Eq. ( p^ . 
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